Rashba split surface states in BiTeBr 
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Within density functional theory, we study bulk band structure and surface states of BiTeBr. 
We consider both ordered and disordered phases which differ in atomic order in the Te-Br sub- 
lattice. On the basis of relativistic ab-initio calculations, we show that the ordered BiTeBr is 
energetically preferable as compared with the disordered one. We demonstrate that both Te- and 
Br-terminated surfaces of the ordered BiTeBr hold surface states with a giant spin-orbit splitting. 
The Te-terminated surface-state spin splitting has the Rashba-type behavior with the coupling pa- 
rameter ur ~ 2 eVA. 

PACS numbers: 73.20.-r, 75.70.Tj 
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INTRODUCTION 

Nowadays, a controllable manipulation of the elec- 
tronic spin degree-of-freedom without recourse to an ex- 
ternal magnetic field is a process in technological de- 
mand, since it constitutes the basis of functionality of 
spintronics devices [l[ . An obvious candidate for the phe- 
nomenon underlying this process is spin-orbit interaction 
(SOI) that couples spin and momentum of electrons. In 
the case of two-dimensional (2D) geometries (surfaces, 
asymmetric quantum wells, etc.), the SOI may result in a 
spin splitting of electron states, which has a nature of the 
so-called Rashba effect (H. This splitting can be tuned 
by an applied electric field 0-0] , which opens a pathway 
for realizing electric-field spin manipulation. There are 
two key operating characteristics here: the Rashba en- 
ergy, Eft, and the momentum offset of split states, k R , 
which together define the Rashba coupling strength as 
a R = 2E R /k R . 

In the conventional semiconductor structures, where 
the Rashba effect has been revealed for the first time, 
the parameter qr is of order of 10 _1 eVA (see, e.g., 
Refs. 0, H, HI). However, for room-temperature appli- 
cations of spintronics devices, it is crucial to have a R 
as large as possible. As a result, over a long period of 
time the Rashba effect keeps attracting a great inter- 
est, and many systems with a large Rashba spin splitting- 
have been discovered. It was found that the f surface 
state on Au(lll) has a Rashba splitting with a R that is 
about five times larger than that in the semiconductor 
heterostructures (see, e.g., Refs. fiol- 12|). A larger a R 
(half as much as that for Au(lll)) has been reported for 
a surface state at the Bi(lll) surface [T3|. In seeking a 
way to tune spin-orbit-splitting of surface states, it was 



shown that, for instance, in the case of the Au(lll) sur- 
face it can be done with deposition of Ag-atoms [lj, [l5| . 
As was demonstrated in Refs. [HHH, a more effective 
way to modify spin splitting of surface states is a surface 
alloying of heavy elements (Bi, Pb, Sb) on noble-metal 
surfaces. In this case, one arrives at a Rashba-type split 
surface state with ccr, that is about one order of magni- 
tude greater than that in the semiconductor structures. 

To go further in possible tuning of spin-orbit-splitting 
of electron states, quantum-well states evolving by the 
confinement of electrons in ultrathin metal films have 
been considered. In the presence of both the surface and 
the interface with a substrate, a number of impacts on the 
splitting doubles. As was reported, e.g., in Ref. a Bi 
monolayer film on Cu(lll) can provide with spin-orbit- 
split quantum-well states in the unoccupied electronic 
structure, which are characterized by cvr similar to that 
in the surface alloys. However, apart from a large spin- 
orbit splitting, for an efficient spintronics application in 
the way specified above a semiconductor substrate and 
the absence of spin-degenerate carriers in a quite wide 
energy interval are more promising. 

In the case of semiconductor substrate (e.g., an ul- 
trathin Pb films on Si(lll) (24[) quantum well states 
show a Rashba splitting as small as in the semicon- 
ductor structures. A large Rashba spin splitting on a 
semiconductor substrate can be reached, for example, by 
means of a Bi-trimer adlayer on a Si(lll) surface (see 
also Ref. [Hj]), where the splitting has a similar origin 
as in the Bi/Ag(lll) surface alloy and a close value for 
the parameter [26|. Nevertheless, the found spin-split 
2D states cannot be well described by a simple Rashba 
model, where a parabolic dispersion with a positive effec- 
tive mass combines with the spin-splitting that is linear 
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in electron momentum. This motivates an active search 
of new materials and a revision of the already known sys- 
tems with a SOI that under certain conditions can lead 
to a technologically meaningful spin spliting of a free- 
electron-like state at a semiconductor surface. In that 
sense, the reexamining of bismuth tellurohalides, where 
a Rashba-type spin splitting of states has been revealed 
to be caused by intrinsic inversion asymmetry of bulk 
crystal potential [U S|, can be considered as a great 
advance made recently in the search. Actually, it was 
shown that Te-terminated surfaces of BiTcCl and BiTcI 
possess a giant s pin- orbit splitting of a free-electron-likc 
surface state 291 32|. 

The mentioned bismuth tellurohalides have hexago- 
nal crystal structures 33| and are characterized by ionic 
bonding with large charge transfer from Bi to halide- and 
Te-atomic layers. The crystal structure is built up of al- 
ternating hexagonal layers Te-Bi-I(Cl) stacked along the 
hexagonal axis. Besides, each three layers Tc-Bi-I(Cl) 
form a three-layer (TL) block, and the distance between 
the blocks is about one and a half times greater than 
the interlayer distances within the Te-Bi-I(Cl) TL struc- 
ture. Such a three-layered structure breaks the inver- 
sion symmetry of bulk crystal potential, which leads to 
appearing of the Rashba-type spin-orbit splitting of the 
bulk bands [27|- Due to the layered crystal structure, the 
bismuth tcllurohalide surfaces can be terminated by Te- 
or halide-atom-layer. Both these terminations hold spin- 
split surface states 2§- 32| . These states emerge by split- 
ting off either from the lowest conduction band (for the 
Te-termination) or from the uppermost valence band (for 
the halide-atom-termination) . The splitting off is caused 
by changes in potential (decreasing at the Te-terminated 
surface and increasing at the halide-atom-terminated sur- 
face) within the near-surface layers [29| as compared with 
the bulk region, which is a consequence of strong ionicity. 

In addition to BiTcCl and BiTel, the bismuth- 
tcllurohalide group is known to have one more 
semiconductor — BiTeBr. It was previously reported 
[33l 34 1 that its layered crystal structure is a disor- 



dered centrosymmetric one, where tellurium and bromine 
atoms randomly distributed within two layers adjacent to 
Bi-atomic layer 3^, 34| . As a consequence, both the bulk 
and surface electronic structure of the disordered BiTeBr 
was never addressed before. Recently, bulk electronic 
structure of the ordered BiTeBr has been calculated and 
Rashba-type splitting of some bands has been analyzed 
& 

In the present paper, we examine both disordered and 
ordered phases of BiTeBr. We model the crystal struc- 
ture of the ordered phase as that of BiTel but with Br 
instead of I and the lattice parameters were taken form 
Ref. [HI], at that atomic positions are obtained within a 
structural optimization. On the basis of ab-initio calcula- 
tions we show that the ordered structure is energetically 
preferable. We demonstrate that both the Te- and Br- 



tcrminated surfaces of the ordered BiTeBr hold the spin- 
orbit split surface states emerged by splitting off from the 
bulk conduction or valence band like in other bismuth 
tellurohalides. For the practical use, as in the case of 
BiTeCl and BiTel the Te-terminated surface of BiTeBr 
is more suitable than the halide-atom-terminated one, 
since it holds a surface state that has a free-electron-likc 
dispersion and a large Rashba-type spin splitting. At 
the same time, BiTeBr has an advantage over BiTeCl 
and BiTel. As compared with BiTeCl, the bromide has 
a larger Rashba spin splitting of the Te-terminated sur- 
face state and a wider bulk band gap. In contrast to 
BiTel, the surface state is larger split off from the bulk 
conduction band and more isotropic. 



COMPUTATIONAL METHOD 

The structural optimization and electronic band calcu- 
lations are performed within the density functional for- 
malism as implemented in VASP [35|,[36| . We use the all- 
electron projector augmented wave (PAW) [13, HH basis 
sets with the generalized gradient approximation (GGA) 
of Perdew, Burke, and Ernzerhof (PBE) [H to the ex- 
change correlation (XC) potential. The Hamiltonian con- 
tains the scalar relativistic corrections, and the SOI was 
taken into account by the second variation method |40j [. 

To treat the bulk disordered phase effect, we employ 
two approaches. The first one is a supercell approach 
used within VASP, where 4x4x1 supercell with several 
configurations for randomly distributed Te and Br atoms 
were considered. The second approach is a virtual crystal 
approximation (VCA) as implemented in the ABINIT 
code 4l|, where the configuration averaged potential of 
a gray atom occupying a site in the Te-Br sublattice is 
defined as a mixture VycA = xVb? + (1 — x)V-[ e of Br 
(^Br) and Te (Vr e ) pseudopotentials with x = 0.5. We 
used GGA-PBE Hartwigsen-Goedecker-Hutter (HGH) 
relativistic norm-conserving pseudopotentials taken 
from Ref. Q which include the SOI. 

The surface of the ordered BiTeBr formed under cleav- 
age can have Te-layer or Br-layer termination. To sim- 
ulate semi-infinite BiTeBr(OOOI), using VASP we con- 
sider a 24 atomic layer slab with bromine side (for Te- 
terminated surface) or tellurium side (for Br-terminated 
surface) passivated by hydrogen monolayer. 



CALCULATION RESULTS AND DISCUSSION 



In earlier works [33j, |3J], the Cdl2 hexagonal struc- 



ture for BiTeBr was reported. This structure differs from 
that of BiTel in that Br and Te atoms are statistically 
distributed over I 2 -type sites (Fig. [lja)). According to 
Ref. [33|, the mixed Te/Br layers are located at a dis- 
tance of ±1.806 A from central Bi layer. To study the 
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FIG. 1: Atomic structure of BiTeBr: disordered structure as 
taken from Ref. [13] (a) and optimized ordered structure (b). 
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FIG. 2: (a) Band structure calculated by the VASP along 
high symmetry directions of the Brillouin zone for the BiTeBr 
ordered phase, and (b) magnified view of the bulk elec- 
tronic structure in the vicinity of the A point calculated with 
use both PAW (VASP) and pseudopotential (ABINIT) ap- 
proaches; inset in panel (b) shows the lowest conduction band 
in the vicinity of the A point calculated within the 4x4x1 
supercell (orange line) and within the VCA (blue line). 



geometry of the disordered phase, we constructed a set 
of 4 x 4 x 1 supcrcells with different configurations of 
randomly distributed Te and Br atoms. The simulation 
shows that owing to markedly different charge transfer 
from Bi atom to Br and Te atoms, the Bi layer acquires 
a substantial rippling of ~ 0.7 A. This value is larger 
than the rippling within the Te/Br layer, which amounts 
to ~ 0.4 A. Such a huge corrugation of atomic layers leads 
to a remarkable variation (from 1.3 to 2 A) of the distance 
between Bi atoms of the central rippled layer and atoms 
of adjacent Te/Br layer. As a consequence, the resulting 
atomic structure differs essentially from that predicted 
in Ref. (33|. Next we considered the BiTcI-type ordered 
structure (Fig.QJb)). In this system, interlayer distances 
obtained within the structural optimization are 1.769 A 
between Bi and Te layers and 1.917 A between Bi and 
Br layers. The ordered BiTeBr gains an energy of ~ 180 
meV per formula unit with respect to the disordered su- 
percells. The VCA calculation confirms the preference of 
the ordered phase. 

In Fig. [21a), we show the bulk band structure calcu- 
lated by VASP code for the ordered BiTeBr. As clearly 
seen in the figure, both the conduction band minimum 
(CBM) and the valence band maximum (VBM) demon- 
strate a giant Rashba-type spin splitting in the H-A-L 
plane. This spin splitting is characterized by a slightly 
anisotropic momentum offset kn that in A-H and A-L di- 
rections is of - 0.05 and 0.04 A^ 1 for the VBM and the 
CBM, respectively. The Rashba energy for the VBM is 
approximately twice of that for the CBM (111 meV vs 66 
meV). As a result, it provides noticeably larger spin-orbit 
coupling for the upper valence-band states as compared 
with the lower conduction-band states (see Tabl. 

As seen in the Fig.^b), the pseudopotential ABINIT 
calculations performed for the ordered BiTeBr confirm 
the large spin-orbit splitting in the vicinity of the A point. 
Moreover, we have obtained values for the Rashba pa- 
rameters, which are very close to those found with VASP 
(see an in Tabl. HJ. The bulk band gap evaluated by 
ABINIT is about 100 meV smaller than that obtained 
from the VASP calculations. The spectrum calculated 
within the VCA for the disordered phase shows practi- 
cally the same band gap and demonstrates the expected 
lack of spin-splitting of the bulk bands due to the pres- 
ence of inversion symmetry in the disordered structure. 
Note that the spin-splitting of the bulk bands obtained 
within supercell approach is negligible, and it agrees well 
with the VCA result (Fig. EJb), inset), which indicates 
that chosen 4x4x1 geometry is well suited for describing 
the disordered BiTeBr. 

As was mentioned above, the spin-split surface states of 
BiTel and BiTeCl emerge as a result of splitting off from 
the lowest conduction band (at the Te-terminated sur- 
face) or from the uppermost valence band (at the halide- 
atom-terminated surface). The splitting off is caused 
by the potential change, AV, in the near-surface lay- 
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TABLE I: Rashba coupling parameters an (eVA) for the bulk 
valence and conduction bands in the vicinity of the A point 
in the A-H and A-L directions. The calculated values for the 
bulk band gap, E g , are also presented. 
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FIG. 3: The change of the potential in the near-surface layers 
of the crystal with respect to that in the central, bulk-like 
layers: (a) Te-terminated surface; (b) Br-terminated surface. 
z = corresponds to the topmost atomic layer. 



ers of the crystal with respect to that in central, bulk- 
like layers. Such a potential change is negative at the 
Te-terminated surface and positive at the C1(I) atom- 
terminated surface 0, [32|, at that AV bears a stepwise 
character owing to clearly defined three-layered structure 
of bismuth tellurohalides. A similar behavior of AV oc- 
curs on the surfaces of the ordered BiTeBr, as seen in 
Fig. [3l where the change of the potential within the three 
outermost TLs on both surface terminations is shown. 

The negative AV observed at the Te-terminated sur- 
face of the ordered BiTeBr leads to a downward shift 
of energies of the electron states trapped in the step- 
wise surface potential (Fig. H|a)). These states are pre- 
dominantly localized in the first three TLs. At the Br- 
terminated surface, an upward shift of energies of the 
trapped states is provided by the positive AV (Fig.[4jb)). 
The trapped states are offset in momentum, reflecting 
a large bulk spin-orbit splitting. They appear partially 
overlapping the valence band continuum except the local 
energy gap regions within bulk continuum states where 
the trapped states can be well resolved in ARPES at high 



binding energies. As a net result, for both terminations 
the changes of the electronic structure of BiTeBr under 
the surface formation lead to emergence of the spin-split 
surface states in the bulk band gap (Fig. [5]). 

At the Te-terminated surface (Fig.[SJa)), the spin-orbit 
split surface state localized in the topmost TL replicates 
the conduction band edge. The degeneracy point of the 
surface spin-split state is 150 meV lower than the CBM. 
Within the energy gap region, the surface-state disper- 
sion demonstrates the free-electron-like parabolic char- 
acter. The spin splitting of the state is characterized by 
a R = 2.109 and 2.007 eVA in the f — K and f — M 
directions, respectively. 

The parabolic character of the surface state provides 
circular shape of the constant energy contours (CEC) for 
inner and outer branches of the spin-split surface state in 
the bulk band gap region. As one can see in Fig.[6ta), in 
approaching the bulk conduction band the CEC for the 
outer branch acquires the hexagonal deformation that is 
already visible at 100 meV above the degeneracy point. 
The surface-state spin structure demonstrates counter- 
clockwise and clockwise in-plane helicity for the inner 
and outer branches, respectively, with a small S z spin 
component for both of them (Fig. Elja)). Owing to sym- 
metry constrains, the expectation value of the S y and S z 
spin components vanishes along F — M , and they have 
maximal values along T — K at any chosen energy. In 
turn, S x is zero along T — K and reaches maximal values 
along T — M direction. In Fig.[6][b), we show the absolute 
value of the cartesian spin components as functions of k\ \ 
for the inner and outer branches of the spin-split surface 
state. As one can see, for a small the \S Z \ component 
is negligibly small, and, thus, the surface state is com- 
pletely in-plane spin polarized. This component starts 
rising at % > 0.1 A -1 , i.e. at energy of 50 meV above 
the degeneracy point, which leads to a decrease of the in- 
plane spin components under approaching the bulk con- 
duction states. Thus, owing to (i) the parabolic energy 
dispersion, (ii) outermost TL localization, and (iii) the in- 
plane helical spin structure within the band gap energy 
region, the surface state on the Te-terminated surface of 
BiTcBr(OOOl) can be described as the Rashba-split sur- 
face state with Rashba coupling parameter of ~ 2 eVA. 

At the Br-terminated surface (Fig. 0(b)), the spin-orbit 
split surface state localized in the topmost TL lies en- 
tirely in the band gap right above the upper valence band 
and follows its edge. As one can see, the states localized 
in the second TL spread already along the valence-band 
edge in T — K and T — M directions but they degenerate 
with bulk states in the close vicinity of T. Furthermore, 
the second outermost-TL localized state arises at f in the 
valley of the valence band. The appearance of the sec- 
ond pair of the spin-split states in the gap and the second 
state at F is explained by the fact that the magnitude of 
the AV is larger than that on the Te-terminated surface 
(sec Fig. [3]). Such a AV provides a larger splitting off 
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FIG. 4: Electronic structure of a BiTeBr(OOOl) slab: (a) Te-terminated surface; (b) Br-terminated surface. The red, pink, and 
orange circles denote weights of the states localized in the 1-st, 2-nd, and 3-rd TLs of the surface under consideration; light 
gray circles mark the states localized on H-terminated side of the slab. The projected bulk band structure is shown in green. 



from the valence band edge. 

In general features, the band-gap-lying outermost 
TL-localized spin-orbit split surface state at the Br- 
tcrminated surface resembles those found at the halide- 
atom-terminated surface of BiTcCl and BiTel [2^, |32| . 
This state demonstrates noticeable anisotropy of the en- 
ergy dispersion with respect to k\\, which results in more 
complex shape of the CECs both above and below the de- 
generacy point (see Fig. [7|). Such an anisotropic disper- 
sion is accompanied by an appreciably out-of-plane spin 
polarization and entangled spin structure of the surface 
state, especially below the degeneracy point (Fig. [7]). 

Formally, the spin-splitting of the Br-terminated sur- 
face state is characterized by kn = 0.079 and 0.077 A -1 
in f — K and f — M direction, respectively, and by Eft 
equal to 130.2 meV (f - K) and 125.6 meV (f - M). 
These characteristics yield or equal to 3.29 and 3.26 eVA 
for T — K and T — M directions, respectively. However, 
this spin-split surface state can not be identified as the 



Rashba state owing to its dispersion and entangled spin 
structure. 



CONCLUSIONS 

Thus, we have investigated the atomic and electronic 
structure of BiTcBr. The total energy calculations of 
the ordered and disordered phases of BiTcBr have shown 
that the ordered structure is energetically preferable. We 
have found that the surfaces of the ordered BiTeBr hold 
surface states which demonstrate a giant spin-orbit spin 
splitting. These states emerge as a result of splitting off 
from the bulk conduction or valence band, owing to the 
potential bending within the near-surface layers, like in 
other bismuth tellurohalides, BiTeCl and BiTel, studied 
earlier. The spin-split surface state at the Te-terminated 
surface, owing to its parabolic energy dispersion, out- 
ermost TL localization, and in-plane helical spin struc- 
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FIG. 5: Magnified view of electronic structure of Te-terminated (a) and Br-terminated (b) BiTeBr(OOOl) surface in the vicinity 
of ? [colors correspond to those marked in Fig. [4]. 



ture preserved within the whole band-gap energy region, 
can be described as a Rashba-split surface state with the 
Rashba coupling parameter of ~ 2 cVA which is in 
a good agreement with recently reported experimental 
value of 2.0(7) eVAjHj]. The Rashba-split state on Tc- 
tcrminatcd BiTcBr has advantages over the other bis- 
muth tcllurohalidcs, which consist in the larger Rashba 
splitting and wider band gap as compared to BiTeCl and 
in the larger splitting off from the bulk conduction band 
with more isotropic energy dispersion within the band 
gap region in comparison with BiTcI. 



The band gap surface state at the Br-terminated sur- 
face owing to its k\ \ anisotropy and entangled spin struc- 
ture can not be identified as the Rashba-split state and 
thus has less appeal than the spin split surface state at 
the Tc-terminated surface. 
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nents S x , S y , and S z at energy of 100 meV (see horizontal 
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